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ABSTRACT 

We present a new technique for directly measuring the distances to superluminal radio sources. 
By comparing the observed proper motions of components in a parsec scale radio jet to their 
measured Doppler factors, we can deduce the distance to the radio source independent of the 
standard rungs in the cosmological distance ladder. This technique requires that the jet angle to 
the line of sight and the ratio of pattern to flow velocities are sufficiently constrained. We evaluate 
a number of possibilities for constraining these parameters and demonstrate the technique on a 
well defined component in the parsec scale jet of the quasar 3C279 (z = 0.536). We find an 
angular size distance to 3C 279 of greater than 1.8^q 3T7 1 ' 8 Gpc, where r\ is the ratio of the energy 
density in the magnetic field to the energy density in the radiating particles in that jet component. 
For an Einstein-de Sitter Universe, this measurement would constrain the Hubble constant to 
be H Q < 65t7 _1/8 km/s/Mpc at the two sigma level. Similar measurements on higher redshift 
sources may help discriminate between cosmological models. 

Subject headings: distance scale - galaxies: active - galaxies: distances and redshifts - galaxies: kine- 
matics and dynamics - quasars: individual (3C 279) - radiation mechanisms: non-thermal 



1. Introduction 

To answer fundamental questions about the ge- 
ometry of the Universe and the distribution of 
matter on the largest scales we must be able to 
measure the distance to objects at large redshift. 
Type la supernovae have begun to fill this role and 
recent results suggest that the overall expansion 
of the Universe may be accelerating (Riess et al. 
f998; Perlmutter et al. 1999). The great bright- 
ness of Type la supernovae, which can out-shine 
their parent galaxies, allow them to be observed 
up to redshifts of z = 1. It is important to have 
cosmological probes at even higher redshifts where 
we are less sensitive to local perturbations (Tipler 
1999) and the predictions of different cosmologi- 
cal models become more distinct (Carroll, Press 
& Turner 1992). 
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Quasars are regularly observed at very large 
redshifts and provide an alternative to gravita- 
tional lenses for high redshift cosmology studies. 
Optically, it is difficult to define a standard candle 
for quasars because of the wide range of intrinsic 
luminosities between objects and the high degree 
of variability of individual objects. Similar difficul- 
ties exist with using the radio luminosity for the 
~ 10% of objects that are radio loud. Global Very 
Long Baseline Interferometry (VLBI) arrays offer 
an alternative to defining a standard candle. The 
ability to resolve compact structures and observe 
proper motions in the radio jets of many sources 
suggests the possibility of defining a standard rod 
from which angular size distances may be derived. 

Kellermann (1993) defined a standard rod for 
compact radio sources as the distance between the 
core and the most distant jet component whose 
peak brightness exceeds 2% of the core brightness. 
He used data from 82 compact sources (out to 
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z w 3) and found an observed relation between an- 
gular size and redshift which is consistent with an 
Einstein-de Sitter Universe (deceleration parame- 
ter, q = 0.5). 

Lynden-Bell (1977) and Lynden-Bell and Liller 
(1978) combined early proper motion measure- 
ments with a light-echo model, e.g. (Couderc 
1939), in the first attempts to use proper mo- 
tions observed in radio jets as cosmological dis- 
tance indicators. Several authors have examined 
the statistics of proper motions of patterns in ra- 
dio jets as a cosmological probe. Yahil (1979) pro- 
posed the idea of a proper motion-redshift diagram 
to measure Hubble's constant and the deceleration 
parameter. Cohen et al. (1988) used measured 
proper motion data and a simple beaming model 
to show that the upper envelope of the proper mo- 
tion data decreased with redshift in a manner con- 
sistent with a Friedmann cosmology and incon- 
sistent with several alternatives. More recently, 
Vermeulen and Cohen (1994) have refined these 
ideas and used a much larger sample to explore 
the possibility of measuring cosmological param- 
eters with proper motion data. They found they 
could usefully constrain cosmological parameters 
and simultaneously learn about the distribution 
of jet parameters. 

Our approach differs from these techniques in 
that we use VLBI observations to directly mea- 
sure the distance to individual superluminal ra- 
dio sources. We compare the proper motions of 
individual components in a parsec-scale radio jet 
with measurements of their Doppler factors. In 
addition, we must constrain the angle the compo- 
nent motion makes with the line of sight and sepa- 
rate the pattern speed (observed in proper motion 
measurements) from the flow speed (observed in 
Doppler factor measurements). In §2 we present a 
number of possibilities for constraining these pa- 
rameters, and in §3 we evaluate the technique with 
an example, using a well defined component in 
the VLBI jet of the quasar 3C 279. Section 4 dis- 
cusses the measurement errors and systematic un- 
certainties associated with this technique. Section 
5 explores the application of such measurements to 
cosmological questions. Our conclusions are pre- 
sented in §6. 



2. Theory 

The observed proper motion, ^, of a pattern 
(component) in a parsec-scale radio jet depends on 
the intrinsic pattern speed, f3 pi the angle between 
the jet axis and the line to the observer, 9, and 
the angular size distance to the radio source, Da- 



j3 p sin I 
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where (3 a is the apparent transverse velocity of the 
component in units of the speed of light. 

We define a Doppler factor for the pattern, 5 P , 

by 



PI 
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The aberrated angle, 8' , is the angle between 
the jet axis and the line to the observer in the 
frame that moves along the jet with the pattern 
speed. It is given by 



cos 9' 
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and 



sin 6' 



dp sin 9. 
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We can use these relations to derive an expres- 
sion for the angular size distance: 
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Equation 5 is plotted for various values of 5 P in 
figure 1. In this equation /1 is in natural units 
of rad/sec; however, in figure 1 /x is in units of 
mas/year. 

The proper motion fi is a directly observed 
quantity. The Doppler factor of the pattern, 5 p , is 
indirectly observed through the use of Synchrotron 
Self-Compton (SSC) or equipartition arguments 
(see §2.2 and the appendix) which measure a prod- 
uct of Doppler factors: 



Sssc = S P (3'f)- 
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and 
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where a is the spectral index (S oc v~ a ) for op- 
tically thin synchrotron radiation and 8'f is the 
Doppler factor of the fluid frame relative to the 
pattern frame: 



f 
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where 6' is the angle of observation in the pat- 
tern frame and is the speed of the fluid in the 
pattern frame. 

For comparison to proper motion measure- 
ments, we are interested in 8 p . To determine 8 P 
from equipartition or SSC arguments we must 
measure or usefully constrain 6f which typically 
requires measuring the pattern versus flow speed; 
however, there is a useful constraint we can place 
on 8'f if we know 9'. For a given 6', the maximum 
in S'f is when [3'f = cos 0' , so 



S' f < 1/y/l -(cOS0') 2 



(9) 



2.1. Additional Constraints 

In the sections that follow, we explore a number 
of constraints that allow us to turn measurements 
of proper motions and Doppler factors into direct 
distance measurements. Some of these techniques 
allow measurement or constraint of S'f which is 
important for accurate determination of S p from 
Doppler factor measurements. 

2.1.1. Bent Jets 

For a given pattern speed, f3 p , the observed 
proper motion is maximized when cos 9 = (3 p . At 
this critical angle, cos0' = and equation 5 re- 
duces to: 
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Sources inside the critical angle will give an up- 
per limit on Da, and sources outside the critical 
angle will give a lower limit. However, these limits 
will only be useful for sources near the critical an- 
gle, and for these sources, we need an additional 



constraint to determine their orientation relative 
to the critical angle. 

Jets which bend on VLBI scales give a unique 
opportunity for observing a source at or near its 
critical angle. As a component on a curved trajec- 
tory passes through the critical angle, a number of 
observable effects occur: the proper motion of the 
component maximizes, the orthogonal component 
of the magnetic field (projected in the plane of the 
sky) maximizes creating a maximum or minimum 
in the observed linear polarization, and thin fea- 
tures, such as shocks, minimize in observed aspect 
ratio. 

We note that if the flow and pattern have differ- 
ent speeds they will also have different critical an- 
gles. The maximization of the proper motion oc- 
curs at the pattern's critical angle. Any maximum 
or minimum in the observed linear polarization de- 
pends on the flow's critical angle. Observing both 
critical angles provides a method of resolving the 
difference between the flow and pattern speeds. If 
only the critical angle of the pattern is observed, 
a useful constraint is that the Doppler factor for 
the flow as observed from the pattern frame, S'f, 
has a maximum value of 1 (see equation 9). 

2.1.2. Aspect Ratio 

Several authors have made use of a sharp (nar- 
row) feature in a radio jet to measure or con- 
strain the jet angle to the line of sight (Eich- 
ler & Smith 1983; Biretta, Owen & Hardee 1983; 
Biretta, Owen & Cornwell 1989; Unwin & Wehrle 
1992). A sharp feature, assumed to be oriented 
perpendicular to the jet direction, is a sign that 
the pattern is moving at close to its critical angle. 
The shape an observer sees for a component mov- 
ing in a radio jet is governed by the aberration 
between the pattern frame and observer frame. 

The observed aspect ratio of a component is 
the ratio of its extent along the jet to its extent 
transverse to the jet , £ = size^/size± (see fig- 
ure 2). The ratio £ constrains cos 9' if the pattern 
is assumed to be axially symmetric and oriented 
perpendicular to the direction of motion. Under 
these circumstances, £ > |cos#'| and we obtain 
potentially useful limits on the angular size dis- 
tance (from equation 5, corresponding to positive 
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and negative signs for cosO 1 ): 
and 

c fJs 2 P +e-i+s P c\ 

C ^KI^)( V J <12 » 

Because the SSC technique produces only a 
lower limit on the Doppler factor (Marscher 1987), 
only the lower limit on Da will be applicable 
when we determine 5 P using that technique. The 
cquipartition assumption will also produce a lower 
limit on the Doppler factor if we only have an up- 
per limit on the frequency of the self-absorption 
turnover in the synchrotron spectrum (Readhead 
1994). 

It is important to note that the relations devel- 
oped here assume that the pattern is perpendicu- 
lar to the jet direction and not oblique. Oblique- 
ness in the plane of the sky is observable from the 
orientation of the feature and perhaps its linear 
polarization. Assuming it is not very large, this 
kind of obliqueness can be corrected for; however, 
if a component can be oblique in the plane of the 
sky it may also be oblique in the plane of obser- 
vation. Obliqueness in the plane of observation 
is indistinguishable from effects of aberration for 
determining the observed component dimensions 
and may cause the relation ( > | cos 9' | to be vio- 
lated. Uncertainty in the degree of obliqueness of 
a given component is equivalent to an added un- 
certainty in the measurement of £ (Birctta, Owen 
& Cornwell 1989; Unwin & Wehrle 1992). 

2.1.3. Linear Polarization 

The linear polarization of a pattern in a ra- 
dio jet provides a measure of the magnetic field 
order. For a tangled magnetic field which has 
been compressed (due perhaps to a propagating 
shock), the degree of linear polarization observed 
depends both upon the degree of compression and 
the viewing angle (Laing 1980; Hughes, Aller & 
Aller 1985). The highest degrees of parallel lin- 
ear polarization (for a given compression) will be 
observed when the flow is moving at or near its 
critical angle. The degree of compression can be 



related to the speeds of the flow (upstream and 
downstream) relative to the propagating shock, 
e.g. (Cawthorne & Wardle 1988; Hughes, Aller 
& Aller 1989). 

Wardle et al. (1994) work out a complete model 
for deducing the jet angle to the line of sight and 
pattern versus flow speeds from detailed VLBI po- 
larization data. They consider the general case of 
a compression in a jet with a tangled field plus a 
component of ordered field along the axis of com- 
pression. By measuring the degree of linear po- 
larization and total intensity in both the shocked 
and un-shocked regions in the jet of 3C 345 (z = 
0.595), they were able to constrain the flow speed 
relative to the shock and the inclination of the jet 
to the line of sight. 

To illustrate the use of the constraints avail- 
able from linear polarization observations, we will 
start with the results of Wardle et al. (1994) for 
jet component C3 of 3C345. The reader is re- 
ferred to that paper for details. For their 1984.2 
epoch, they find nominal values of 6 = 2.5°/Da, 
(3d = 0.5, and (3 a = 11. IDa (m = 0.44 mas/year). 
The flow speed of the shocked fluid towards the 
core in the frame of the shock- front is (3d- In the 
notation of this paper, (3'j = —(3d- Da is measured 
in Gpc for the values given above. 

We use equations 3 and 7 to calculate cos 8' = 
0.5 and 5'f — 0.7 which are essentially independent 
of distance. If we had observations of the syn- 
chrotron self-absorption turnover for component 
C3 in 3C345 at this epoch, we could measure its 
total Doppler factor using equipartition or SSC ar- 
guments. We could then use the measurement of 
6'f to determine the Doppler factor of the pattern, 
5 p . With dp, n, and cos#' determined, equation 5 
would allow calculation of the distance to 3C345. 

2.1.4- Jet/ Counter- Jet Ratio 

Another potentially useful constraint is the ob- 
served jet/counter-jet brightness ratio, i?, e.g. 
(Unwin & Wehrle 1992). 

where n = 3 for discrete components and n = 2 
for continuous jet emission. The reduction of the 
equation to include (3 a is only valid if the pattern 
and flow speeds are the same. If they are the same 
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or if we know the relationship between them, mea- 
surement of R allows us to directly compare the 
Doppler factor to observed proper motion and de- 
duce the distance to the source. In the event that 
the pattern and flow speeds are the same, the an- 
gular size distance is given by 



D A = 



R— 



8 2 
v 



(14) 



An attractive feature of this approach is that 
the final answer depends weakly on the measure- 
ment of R. For highly beamed sources, how- 
ever, R is huge and even high quality VLBI mea- 
surements cannot usefully constrain it. For less 
beamed sources, measuring or constraining R is 
more promising. A major drawback to using R to 
connect fi and S p is that R is a global property of 
the jets rather than of an individual component. 
For R to be useful, the flow speed and angle to 
the line of sight need to be the same for the jet 
and counter-jet and constant over the region for 
which R is measured. In addition, there cannot 
be significant excess absorption of emission from 
the counter-jet, c.f. (Krichbaum et al. 1998). 

2.2. Doppler Factors 

The synchrotron spectral turnover provides a 
kind of natural (but broad) spectral line for ho- 
mogeneous synchrotron sources. By carefully mea- 
suring the spectrum of a component and its angu- 
lar size we can use limits on its SSC x-ray flux, e.g. 
(Marscher 1987) or an assumption of equipartition 
between the field and particle energies (Readhead 
1994) to determine the Doppler factor. In the ap- 
pendix these formulae are presented for arbitrary 
homogeneous geometry and for the specific case of 
a spherical geometry. 

We have chosen to use a spherical component 
geometry for the calculations presented in this pa- 
per. Without detailed knowledge of the true geom- 
etry a spherical geometry is well-suited for calcu- 
lation because it computes the angular area pre- 
sented to the observer reasonably, gives a sensi- 
ble line of sight through the component, and pro- 
vides naturally for a range of optical depths across 
the component. Marscher (1987) suggests using 
9 d ~ 1.8y/e Ga 9 Gb to convert Gaussian FWHM di- 
mensions (measured in model-fitting) to spherical 
diameters, and we adopt his approximation. (See 



appendix B for discussion of the effect of assumed 
model geometry on measured Doppler factor.) 

3. Example: 3C 279 

As an example of these ideas, we use a well de- 
fined component in the milli-arcsecond jet of the 
well known blazar 3C279 at z = 0.536. We ob- 
served 3C 279 with the Very Long Baseline Array 3 
(VLBA) for six epochs at 15 and 22 GHz during 
1996 and at four frequencies (5.0, 8.4, 15.4, and 
22.2 GHz) during December of 1997. These obser- 
vations were all calibrated using standard tech- 
niques, e.g. (Cotton 1993; Roberts, Wardle & 
Brown 1994), using the National Radio Astron- 
omy Observatory's Astronomical Image Process- 
ing System (AIPS). Model-fitting was performed 
in the (u,v)-plane with the Caltech VLBI program, 
DIFMAP. 

3.1. Proper Motion 

Figure 3 shows the structure of the inner jet of 
3C279 at 22 GHz. Table 1 gives detailed com- 
ponent data from model-fitting the inner jet of 
3C279 for the 1997.94 epoch. Kl is a well de- 
fined, strong component which has persisted for 
years, e.g. (Unwin et al. 1998). The component is 
located approximately 3 milli-arcseconds from the 
core at a position angle of —115°. Over the course 
of our observations we observe this component to 
move radially from the core with a proper mo- 
tion of fi — 0.24±0.01 mas/year (see Figure 4). It 
maintains a structural position angle of — 114° ±1° 
over the course of our observations. 

3.2. Synchrotron Self- Absorption 
Turnover 

We fit the total intensity of Kl at all four fre- 
quencies in 1997.94 and have measured its spec- 
tral turnover. Figure 5 displays the fit of a syn- 
chrotron self-absorption spectrum to the data as- 
suming a slab geometry. (Fitting the spectral 
shape of a homogeneous sphere gives a nearly iden- 
tical result but a slightly smaller error range on 
the parameters.) The spectral turnover is located 



3 The VLBA is part of the National Radio Astronomy Obser- 
vatory, which is a facility of the National Science Founda- 
tion operated under cooperative agreement by Associated 
Universities, Inc. 
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at v peak = 6.02 GHz (+0.33, -0.49) with a flux, 
Speak = 4.40 Jy (+0.19,-0.07), and a spectral 
index, a = 0.52 ± 0.05. The errors in the fit 
are approximately 1 a errors found by a Monte 
Carlo simulation. The simulation created and fit 
1000 fictional data sets using the measured data 
and assuming the measurements are Gaussian dis- 
tributed with 1 a deviations given by the measured 
error bars. The error bars on the fluxes were esti- 
mated by varying parameters in the model-fits us- 
ing Jim Lovell's Difwrap program, an interactive 
shell for the Caltcch VLBI program, DIFMAP. A 
number of factors were used to gauge the size of 
the error bars including shape of the chi-square 
minimum, noise on the residual map, and direct 
comparison of model and data in the (u,v)-plane. 
We have no direct way of knowing if the errors 
estimated for the fluxes are genuinely 1 a errors; 
however, the spectral fit has 1 degree of freedom 
(4 data points and 3 parameters) so the x 2 of the 
spectral fit should be near unity if the errors on 
the data are 1 a. The measured x 2 of the spectral 
fit is 1.0. 

The chief uncertainty in the measured fluxes of 
Kl is due the presence of its poorly defined "tail", 
fit as component K2. While Kl is fit robustly by a 
sharp Gaussian component, K2 is more difficult to 
fit. This becomes more of a problem at the lower 
frequencies (especially 5 GHz) where Kl is not as 
well resolved. 

To check the spectral fit for Kl, we examined 
the spectrum of observed fractional linear polar- 
ization. Figure 6 displays the observed fractional 
polarization plotted together with a theoretical 
curve produced by numerical simulation. The sim- 
ulation is of a homogeneous slab with the same 
total intensity spectrum as fit to Kl. The simu- 
lation solved the full equations of polarized trans- 
fer, e.g. (Jones & O'Dell 1977), for a completely 
tangled magnetic field plus a small ordered com- 
ponent. The magnitude of the ordered component 
was scaled to give the observed fractional polar- 
ization at 22 GHz. To simplify the simulation, no 
internal Faraday rotation was allowed. It is clear 
that the spectrum of the fractional polarization 
is completely consistent with total intensity spec- 
trum fit to Kl. 



3.3. Doppler Factor from Equipartition 

We now use the measured angular size of Kl 
at 22 GHz and the fit to the spectral turnover to 
deduce an equipartition Doppler factor (derived in 
the appendix). The measured FWHM angular size 
at 22 GHz is 9 Ga x9 Gb = 0.46(±0.02)x0.20(±0.01) 
mas x mas. The error bars were estimated by 
varying parameters in the model-fit and by com- 
parison to the measurements at 15 GHz and 8 
GHz. We use the Doppler factor formulation for 
a homogeneous sphere (equation A6) with 9d ~ 
1.8y/6 G J Gb . So for Kl, 6 d = 0.55 ± 0.02 mas and 
we obtain an equipartition Doppler factor for the 
pattern of 

s P = i9.i±E:i V7 (ty" 6 ' 7 (15) 

where r\ = Ub/U tp is the equipartition factor 
(Ub = magnetic field energy density; U rp = energy 
density in the radiating particles) , Da is the angu- 
lar size distance in Gpc, and S'j is the Doppler fac- 
tor of the flow as viewed by an observer co-moving 
with the pattern. To do this computation, we have 
assumed energy spectrum limits of — 10 and 
72 = 1 x 10 6 . For a ~ 0.5, the dependence on 
these limits is approximately [ln(7 2 /7i)] 1 ^ 7 . 

It is interesting to compare this Doppler factor 
for Kl to a Doppler factor measured for the com- 
ponent K4. The spectral fit for K4 is given in Fig- 
ure 7, we find v peak = 12.59 GHz (+0.70, -0.41), 
S peak = 10.86 Jy (+0.40,-0.37), and a = 0.70 
(+0.17,-0.16) with x 2 = 0.8. We have only an 
upper limit on the angular size of the component 
transverse to the jet direction, so we can only 
use the limit 9 d < 0.42 ± 0.01. We calculate an 
equipartition pattern Doppler factor, 

S P > 16.7±S;S (^-) 1/7 4 {5> } y^ (16) 

for the component K4. The dependence on the 
energy spectrum limits for a ~ 0.7 is approxi- 
mately 7^ ' 4 / 7 ' 4 . Ghisellini et al. (1993) report 
an SSC Doppler factor for the core of 3C 279 of 
Sssc > 18.0. 

3.4. Aspect Ratio 

Kl is a narrow component oriented perpendic- 
ular to its position angle. At 15 and 22 GHz the 
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component is well resolved in both directions, at 
8 GHz it is less well resolved, and at 5 GHz the 
component is unresolved along the direction of the 
jet. Assuming that the component is not oblique, 
the observed aspect ratio is ( = 0.43 ± 0.03. 

Component Kl shows no sign of significant 
obliqueness. At 22 GHz we measure its orien- 
tation to be 7° ± 11° from perpendicular to its 
long term structural position angle. The high fre- 
quency linear polarization of Kl is aligned with its 
long term structural position angle to —3° ± 6°. 
(The overall calibration of the polarization posi- 
tion angle for our VLBA observations in Dec. 1997 
was from simultaneous VLA observations of the 
compact source OJ287.) Using these estimates on 
the obliqueness in the plane of the sky as a guide, 
we estimate an uncertainty in the obliqueness in 
the plane of observation of ±5°. This uncertainty 
in the degree of obliqueness translates to an ad- 
ditional uncertainty in the measured aspect ratio, 
roughly ( = 0.43 ± 0.08. 

3.5. Measuring the Distance 

Using our measurement of the observed proper 
motion, pattern Doppler factor, and aspect ratio 
of component Kl, equation 11 gives the following 
limit on the angular size distance to 3C279: 

D A >l.St° iv 1/S Gpc (17) 

which depends only on the equipartition factor, 
i]. Note that we have used equation 9 and our 
measurement of ( (> | cos#'|) to limit the Doppler 
factor of the flow relative to the pattern frame to 
S'j < 1.1 ± 0.1. (Because we have only an upper 
limit on S'j, the upper limit on the distance (equa- 
tion 12) is undetermined.) 

4. Discussion 

Equation 17 provides only a lower limit on the 
angular size distance to 3C279. The result is a 
limit because we could only make use of the as- 
pect ratio constraint. In general this technique 
can provide direct measurements (not just limits) 
for sources where some of the other constraints in 
Section 2.1 can be applied successfully. 

The 1 a errors on this limit are +28% and 
— 17%. These errors are dominated by the un- 
certainty in the spectral turnover measurement 
of component Kl. For Kl, we have only one 



spectral point on the optically thick side of the 
turnover and this point is poorly constrained 
due reduced resolution at 5 GHz. The spectral 
turnover for component K4 is better determined 
giving a Doppler factor with 1 a errors of +22% 
and —13%. With better frequency coverage (per- 
haps by using widely separated IF channels near 
the spectral turnover) and better angular resolu- 
tion at lower frequencies (through the use of space 
VLBI), we believe we can eventually reduce the 
1 a measurement errors on Doppler factors from 
equipartition and SSC techniques to ~ 10 — 15%. 

For our 3C 279 distance limit, the equipartition 
factor, 77, is the most significant unknown quan- 
tity. Even though 77 enters to only a small fac- 
tor, it is poorly constrained. Readhead (1994), 
when proposing the technique, argued that sources 
should be near equipartition (77 = 1 for electron- 
positron jets) and suggested an error of ~ 13% 
in the Doppler factor for typical departures from 
equipartition. Singal (1986) calculated the dia- 
magnctic effect of spiraling electrons in a magnetic 
field. He found that the energy density of the elec- 
trons could not exceed 6 times the energy of the 
magnetic field and still maintain synchrotron ra- 
diation. Bobo, Ghisellini & Trussoni (1992) re- 
peated Singal's calculation and included a surface 
current term. They found that the energy density 
of the electrons could not exceed the energy den- 
sity of the applied field by more than a factor of 3 
although they note that the energy density of the 
effective field in the region could be much smaller 
than the applied field. 

Another issue important to calculating pre- 
cise values for equipartition Doppler factors is 
the cutoffs in the power law particle energy spec- 
trum. The computed Doppler factors depend only 
weakly on the assumed value of the cutoffs, but 
different, reasonable assumptions for the cutoffs 
may lead to a <~ 5 — 10% uncertainty in the com- 
puted Doppler factor. Energy spectrum cutoffs 
can be measured or constrained, however. In War- 
die et al. (1998) we used VLBI circular polariza- 
tion observations to show the lower energy cutoff 
in 3C 279 was 71 < 20. The uncertainty in the low 
energy cutoff dominates when a > 0.5. 

One way around the uncertainties in assum- 
ing equipartition and energy spectrum cutoffs is 
to calculate Doppler factors using measured x-ray 
fluxes. A drawback to this approach is that the 
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SSC Doppler factor calculated for a given compo- 
nent is only a lower limit because, with the ca- 
pabilities of current instruments, observed x-ray 
fluxes include contributions from all parts of the 
parsec-scale source. Such a limit is unlikely to be 
useful for jet components like Kl in 3C 279 which 
contribute a very small fraction of the total x-ray 
flux of the source. Assuming rj = 1 and invert- 
ing equation A9 to solve for the x-ray flux of Kl 
yields ~ 0.01/zJy of 2 KeV x-rays. This flux is 
less than 1% of the total 2 KeV x-ray flux re- 
ported by Wehrle et al. (1998) for 3C279 in its 
quiescent state in January of 1996. It is interest- 
ing to note that the capabilities of the proposed 
MAXIM program (http://maxim.gsfc.nasa.gov/) 
would make x-ray observation of individual jet 
components possible. 

The final area of systematic uncertainty is the 
assumed geometry for the pattern. We found that 
we obtained essentially identical spectral fits when 
using the functional form for a uniform slab as for 
a uniform sphere, so we can safely say that the as- 
sumed geometry has little affect on the spectral fit. 
In appendix B we explore the remaining depen- 
dence of derived Doppler factors on assumed pat- 
tern geometry. The main result is that a spherical 
geometry should be a good approximation for cal- 
culational purposes and tends to produce a lower 
limit on the Doppler factor if the true geometry is 
non-spherical. 

5. Application to Cosmology 

The angular size distance in terms of redshift, 
z, Hubble constant, Hq, matter density, SIm, and 
cosmological constant, = A/(3Hq) is given by, 
e.g. (Carroll, Press & Turner 1992) 

Da = j= 

x s(<v/R / [(i + *') 2 (i + JW) 

Jo 

-z'(2 + z')Q A }- 1 ^dz') (18) 

where if Qm + > 1 then §(x) = sin(x) and k = 
1 — f!ju — Oa, if f^M + f^A < 1 then §(x) = sinh(x) 
and k = 1 — Om — ^A, and if Qm + Cl\ = l then 
S(x) = x and k = 1. 

In figure 8 we plot our lower limit on the dis- 
tance to 3C 279 (assuming i) = 1) against two cos- 
mological models for a flat universe (Am + Oa = 



1). The first case is consistent with recent type la 
supernovae results (Perlmutter et al. 1999; Ricss 
ct al. 1998). The second case is the standard 
Einstein-de Sitter universe. The Hubble constant 
is maintained as a free parameter which scales the 
ordinate of figure 8. 

Our single distance limit cannot distinguish be- 
tween these cosmological models; however, we can 
investigate whether this technique holds promise 
for distinguishing these cases in the future. At 
redshifts larger than 1.0 these cosmological models 
differ by as much as 30 — 40%. Superluminal radio 
sources are regularly observed at these large red- 
shifts. For 3C 279 our lower limit has 1 a measure- 
ment errors on the order of 20—25% and we believe 
it is reasonable obtain distance measures (or lim- 
its) good to ~ 10 — 15% for carefully planned ob- 
servations on well selected objects. Even a handful 
objects (over a range of redshift values) could pro- 
vide strong constraints on cosmological models. 

The systematic uncertainties (discussed in §4) 
in measuring the Doppler factor present the 
largest difficulty here. We can constrain the range 
of allowed values for the equipartition factor, r], 
by comparing distance measurements to sources at 
similar redshifts. To detect any systematic offset 
in 7] from unity, we must have a way of calibrating 
our Doppler factor measurements. Fortunately, 
all of the poorly constrained quantities (77, en- 
ergy spectrum cutoffs, geometric dependence) can 
be grouped as a single multiplicative parameter in 
the equipartition Doppler factor. We can calibrate 
any systematic offset in this parameter by mak- 
ing distance measurements to sources at moderate 
redshifts (where the effects of differing cosmolog- 
ical models are not strong) and comparing the 
result to other distance measurement techniques. 
Such a scheme would introduce a dependence on 
the cosmological distance ladder, so we should 
actively seek other techniques for calibrating sys- 
tematic effects in Doppler factor measurements. 

6. Conclusions 

We have demonstrated a new technique for di- 
rectly measuring the distances to high redshift, su- 
perluminal radio sources. This technique involves 
the comparison of Doppler factor and proper mo- 
tion measurements for individual jet components; 
we must also determine the jet angle to the line 
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of sight and pattern versus flow velocity. We 
have presented several techniques for measuring 
or usefully constraining these parameters. In gen- 
eral these techniques will be applicable only to 
selected sources which have jet components with 
the right characteristics; however, with hundreds 
of currently known supcrluminal sources and new 
jet components emerging frequently from many of 
them, it seems reasonable to assume that we will 
find a significant number of candidates. One inter- 
esting possibility is that some sources could have 
more than one component to which we can apply 
these techniques, giving us multiple, independent 
distance determinations to the same object. 

To begin to usefully constrain cosmological pa- 
rameters we need to obtain high quality distance 
measurements or limits to several sources over a 
range of redshift. We have performed a detailed 
analysis of the measurement error associated with 
our limit on the distance to 3C 279. We found the 
measurement error was ~ 20 — 25% and concluded 
that carefully planned, high-quality observations 
could reduce measurement error to <~ 10 — 15%. 
The systematic error was more difficult to quan- 
tify, although we outlined the sources of system- 
atic error and presented rough estimates. Apply- 
ing this technique to a larger sample of sources will 
be important not only for probing cosmological pa- 
rameters but also for investigating the sources of 
systematic error and how close these sources are 
to equipartition. 
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A. Doppler Factor Formulae 

While the relativistic plasma in radio jets does not contain atoms and molecules whose characteristic 
spectra we could use to directly measure their bulk Doppler factors, they do contain a unique (but broad) 
spectral feature: the synchrotron self-absorption turnover. The location of the turnover due to synchrotron 
self-absorption depends on the magnetic field strength, particle density, and size of the emission region. For 
a volume of homogeneous plasma, we can use the extrapolated optically thin flux at the turnover frequency 
and the optical depth at the turnover frequency to solve for the magnetic field and the particle density, e.g. 
(Marscher 1987). 



B = jy^Cst^S-W (Al) 

and 

where K is the constant in the power law particle density, N^dj = K^^ 2a+1 ^d^. The sign of the spectral 

index, a, for optically thin emission is given by S ex v~ a . The parameter Q — Jf , and £ c = s' c /Da, 

where V is the volume in the pattern frame and s' c is the line of sight through the center of the volume 
along which the optical depth at the turnover, r m , is defined. The extrapolated optically thin flux 4 at the 
turnover frequency, v m , is given by S m . Constants C3 and C4 are tabulated in table 2 for £1 in mas 2 , v m in 
GHz, S m in Jy, Da in Gpc, and £ c in mas. 

Assuming we can translate the source dimensions, V 1 and s' c , into observable quantities (e.g. observed 
angular size), these expressions provide us with two equations and three unknowns, B, K, and the Doppler 
factor. To solve for these quantities, it is necessary to have a third constraint. Two possibilities for a third 
constraint are equipartition (Readhead 1994) and Synchrotron Self-Compton (SSC) x-ray flux (Marscher 
1987). 

A.l. Equipartition 

Equipartition assumes that the energy density of the magnetic field and the energy density of the particles 
are equal. We will parameterize the relationship between the energy density of the magnetic field and the 
energy density in the radiating particles: 

U B = vU rp (A3) 

where Ub — B 2 /Sir, U rp — j mc 2 jN 7 d'y. For an electron-positron jet, r\ = 1 for equipartition. In an 
electron-proton jet, the value of n for equipartition will depend on the details of the particle acceleration 
within the jet. (See §4 for discussion of the value of 77.) 
The equipartition condition is 

B 2 /8n = r 1 mc 2 g(a, 71, j 2 )K, (A4) 

where for a ^ 0.5 

/ \ 1 ( -(2a-l) -(2a-l)\ 

g(a,7i>72) = 2a - 1 V 1 V ) 

4 The factor by which S m over-predicts the observed peak flux, S Q , is tabulated in table 2 for a spherical geometry. 
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and for a = 0.5 



ff(a,7i,72) = In — 
\ 7i 



71 and 72 represent the lower and upper cutoffs for the particle energy distribution. 

Substituting the expressions for B and K into the equipartition condition and solving for the Dopplcr 
factor yields 



8 eq =8 p {5' f )^=F eq r-\l + z) 



ff(a, 71, 72)775 



2a+7 



2a + 6 



where F eq is tabulated in table 2 for SI in mas 2 , v m in GHz, S m in Jy, D A in Gpc, and £ c in mas. 



(A5) 



Homogeneous Sphere Geometry: The spherical geometry assumes that the components emitted from 
the AGN are homogeneous spheres of radiating plasma. The sphere has a radius, R, and we can define 
an angular diameter, 0d — 2R/Da- Therefore ft — (ir/6)9 d and £ c = 64. With these identifications, the 
equipartition Doppler factor is given by 



6 eq = F eq T m 1 (l + z) ( £ 

7T 



»+7 
* + 6 



g(a, 71, 72)775; 



2a+7 



D A (e d y m ) 4a+15 



(A6) 



The equipartition Doppler factor derived by Readhead (1994) has a slightly different functional form than 
our expression. Readhead's expression is derived by comparing observed brightness temperature, T b ', to an 
(rest-frame) equipartition brightness temperature, T eq : S eq — T b '/T eq . For 3C 279 we can turn this expression 
around and calculate T eq for component Kl from our measured S eq and observed brightness temperature. 
We find T eq = 4 x lO 10 ^ (for 77 = 1) which is close to T eq ~5x 10 10 K which Readhead argues is a typical 
upper cutoff for powerful extra-galactic radio sources in their rest-frame. 

A. 2. Synchrotron Self-Compton Emission 

A second possible constraint is synchrotron self-Compton x-ray flux, e.g. (Marscher 1987). The observed 
x-ray flux from the SSC process is given by (adapted from Rybicki & Lightman (1979)) 



(A7) 



where 



A(p) = 2 p+3 



4p+ 11 



(p + 3) 2 (p + 5)(p+l) 



p = 2a + 1 



and 5 is a factor (<~ 1) that accounts for the differences in photon number density (resulting from edge 
effects) within the emitting volume, I" is the average line of sight as seen from the center of the geometry in 
the fluid frame, ff, is the upper cutoff frequency of the synchrotron emission spectrum, and v c is the x-ray 
observation frequency. As defined earlier, S m is the extrapolated, optically thin synchrotron flux at the 
turnover frequency, v m . 

Substituting for K and solving for the Doppler factor gives 



Sssc = S p (6' f )£** = FsscS m (l + z) 



"in In 



-(3a+5) -(2a+2) 
L/rn I m 



Sx{hv c ) a KeV W«+* 



(A8) 
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where Fssc is a constant which is tabulated in table 2 for fl in mas 2 , v m in GHz, S m in Jy, and Sx in A*Jy- 
It is important to note that, in practical application, the SSC Doppler factor is only a lower limit. With 
the capabilities of current instruments, observed x-ray fluxes include contributions from all parts of the 
parsec scale jet, thermal x-ray emission from the accretion disk, and inverse-Compton x-rays which are not 
the result of the SSC process. We therefore obtain only an upper limit on the SSC x-ray flux of a given 
component and thus a lower limit on its Doppler factor. 

Homogeneous Sphere Geometry: O = (7r/6)6>^ and s' c = 2R. Also, for a spherical geometry 5 I" = R 
and S = 3/4 (Gould 1979). Making these substitutions gives the following expression for the SSC Doppler 
factor: 



B. Choice of Model Geometry for Doppler Factor Measurements 

Differences between model geometries show up in the range of optical depths across the component (which 
affects both r m and the ratio S m /S ) and conversion of measured Gaussian FWHM diameters, Oq, to the 
angular dimensions of the assumed geometry. These conversion factors can be estimated for simple geometries 
by matching the second moment of their Fourier transforms which are fit in the (u,v)-plane by Gaussian 
models. For a uniformly bright disk, Od — 1.7#g- For a homogeneous sphere, Od — 1.90G- 

For the purposes of measuring the Doppler factor, the effect of having a wide range of optical depths for a 
homogeneous sphere nearly offsets the larger conversion factor for Gaussian measured component dimensions. 
If we assume that a given component is a sphere, when in reality it is a uniformly bright disk, we will calculate 
a Doppler factor that is about 10% too small using the equipartition formula or about 5% too small using 
the SSC formula. In this scenario there is an additional factor by which assuming a sphere will under-predict 
the Doppler factor. This factor is due to the unknown physical depth of the uniformly bright disk. For the 
equipartition formula, this factor comes in as <~ {Od/i) 1 ^ 1 , where £ is the angular thickness of the disk. For 
the SSC formula, this factor is - [(8/3)(l" /s')E] 1/5 , where I" is the mean line of sight as seen by a photon 
at the center of the disk (in the fluid frame) , s' is the physical depth of the disk (pattern frame) , and 5 <~ 1 
accounts for differences in photon density throughout the disk. 

In general, the patterns we observe are likely to be some compromise between these geometries, perhaps a 
cylindrical disk viewed nearly edge on. Without detailed knowledge of the true geometry a spherical geometry 
is well-suited for calculation because it computes the angular area presented to the observer reasonably, gives 
a sensible line of sight through the component, and provides naturally for a range of optical depths across 
the component. Marscher (1987) suggests using 6d — 1.8-v/^g ^G(, to convert Gaussian FWHM dimensions 
to spherical diameters, and we have adopted his approximation. 



This is strictly true only if the pattern and flow are moving with the same speed. 




(A9) 
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Fig. 1. — A plot of angular size distance in Gpcx^i(l + z) versus cos 9' where 9' is the angle of observation in 
the pattern frame and /i is the observed proper motion in mas/yr. The very large distances near cos 9' = — 1 
are a curious feature of this plot. Patterns with a negative cos#' are viewed from behind and travel at a 
larger angle than the critical angle, /? = cos6* c . A source which has a high Dopplcr factor, 8 P , and a negative 
cos 9' must have a very small critical angle, 9 C . With a very small critical angle, such a source has a very 
large (3 giving it a large apparent speed, f3 a oc Da x fi. 




Fig. 2. — An observer sees the profile of an axially symmetric pattern from the aberrated angle, & . The 
measured aspect ratio, £, is always greater than or equal to cos& . 
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Fig. 3.— Naturally weighted image of the jet of 3C279 at 22 GHz. Epoch 1997.94. The locations of the 
components Kl, K2, K4, and D are marked on the image. Component data are summarized in table 1. 



Proper Motion of Kl 
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CD 
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Fig. 4. — Proper motion of component Kl. Component positions are taken from the 15 GHz model-fits. 
With the data weighted equally, the derived proper motion is = 0.24 ± 0.01 mas/year. Error bars on the 
positions are not plotted because we do not have a good a-priori method for estimating them. From the 
small deviation of the data from the fit, we can say that the errors in the radial position of Kl are ~ 0.02 
which is about 1/20 th of the uniformly weighted beam width along that direction. 
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Spectral Fit for Kl 
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Frequency (GHz) 

Fig. 5. — Synchrotron self-absorption spectrum of component Kl. v pea k = 6.02 GHz (+0.33, —0.49), S pea k = 
4.40 Jy (+0.19, -0.07), and a = 0.52 ± 0.05. 



Polarization Spectrum for Kl 




Frequency (GHz) 

Fig. 6. — Fractional linear polarization of Kl plotted against frequency. The "theoretical curve" is from a 
numerical simulation of a homogeneous slab with a completely tangled magnetic field plus a small ordered 
component. The size of the ordered component was scaled to match the observed fractional polarization at 
22 GHz. The opacity of the slab was fixed to match the fit to the total intensity spectrum. 
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Spectral Fit for K4 




Frequency (GHz) 

Fig. 7. — Synchrotron self-absorption spectrum of component K4. Vpeak = 12.59 GHz (+0.70,-0.41), 
S peak = 10.86 Jy (+0.40, -0.37), and a = 0.70 (+0.17, -0.16). 



Distance vs. Redshift 
W I 1 1 1 1 1 1 , , , 1- 




Redshift 



Fig. 8.— Angular size distance versus redshift for two cosmological models for a flat universe (Om + ^a = !)• 
The fijif = 0.3, fl\ — 0.7 case is consistent with recent SNe la results, and the Qm — 1-0, Oa = 0.0 case 
corresponds to an Einstein-de Sitter universe. The Hubble constant is taken to be Hq = lOO/i km/s/Mpc. 
Our distance limit to 3C279 (assuming r\ = 1) is plotted for various values of h. The points for h = 0.65 
and h = 0.80 are offset slightly in redshift to enhance readability. 
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Table 1: Multi-frequency Component Data for 3C279 for epoch 1997.94. 



Component 


Frequency 




niL 


X 


R 


e 


ma j .axis 


min.axis 


4> 




(GHz) 


(Jy) 


(%) 


(dog) 


(mas) 


(dog) 


(mas) 


(mas) 


(dog) 


D 


4.99 


5.20 


8.9 


-82 














8.42 


7.31 


7.4 


-46 














15.37 


12.08 


10. 6 a 


_4 a 














22.23 


16.27 


8.2 a 


-10 a 






0.10 


< 0.15 


72 


K4 


4.99 


2.51 


3.1 


-28 


0.48 


-121 










8.42 


8.84 


2.9 


41 


0.33 


-126 


0.38 


< 0.50 


61 




15.37 


10.42 


10. 9 a 


84 a 


0.32 


-132 


0.37 


< 0.25 


71 




22.23 


8.96 


10. 2 a 


70 a 


0.35 


-135 


0.36 


< 0.15 


70 


K2 


4.99 


0.62 


17.5 


81 


2.88 


-123 


1.06 


< 0.30 


-17 




8.42 


0.29 


21.2 


88 


2.87 


-130 


2.08 


0.21 


-30 




15.37 


0.23 


17.4 


81 


3.12 


-124 


1.23 


0.27 


7 




22.23 


0.19 


11.8 


66 


3.07 


-125 


1.31 


< 0.07 


2 


Kl 


4.99 


4.26 


4.2 


93 


3.21 


-114 


0.41 b 


< 0.30 


-24 




8.42 


4.17 


6.0 


76 


3.31 


-114 


0.46 b 


0.18 b 


-27 




15.37 


3.11 


5.9 


64 


3.36 


-115 


0.44 


0.19 


-19 




22.23 


2.68 


6.5 


63 


3.37 


-115 


0.46 


0.20 


-17 



Note. — Limits on angular size are estimated to be 1/5 of the uniformly weighted beam width along the corresponding 
dimension. 

"These values were obtained by allowing the position of the core, D, to float in linear polarization and should be used with some 
caution. 

'These angular dimensions appear resolved in the model-fit but are somewhat smaller than the formal limits. 



Table 2: Tabulated Constants. 



a 






C 3 




c 4 


F eq 


Fssc 


T m 


S„ 


./So 





.1 


2 


14 


X 


10 


-a 


1.50 


1.36 


1.02 


0.104 


1 


.04 





.2 


1. 


.39 


X 


10 


-3 


10.2 


2.08 


1.01 


0.204 


1 


.08 





3 


9 


.62 


X 


10 


-4 


66.3 


3.02 


1.00 


0.300 


1 


.12 


0. 


.4 


7. 


.01 


X 


10 


-4 


417 


4.20 


0.987 


0.392 


1 


.15 


0. 


.5 


5. 


.30 


X 


10 


-4 


2.57 x 10+ 3 


5.67 


0.972 


0.480 


1 


.19 





.6 


4 


12 


X 


10 


--1 


1.56 x 10+ 4 


7.44 


0.956 


0.565 


1 


.22 





.7 


3 


.28 


X 


10 


-4 


9.42 x 10+ 4 


9.56 


0.942 


0.648 


1 


.26 





.8 


2. 


.66 


X 


10 


-4 


5.65 x 10+ 5 


12.1 


0.929 


0.727 


1 


.29 


0. 


.9 


2. 


.19 


X 


10 


-4 


3.38 x 10+ 6 


15.0 


0.918 


0.804 


1 


.33 


1. 





1. 


.82 


X 


10 


-4 


2.02 x 10+ 7 


18.3 


0.907 


0.878 


1 


.36 


1. 


.1 


1. 


.54 


X 


10 


-4 


1.21 x 10+ 8 


22.1 


0.899 


0.951 


1 


.39 


1 


.2 


1 


31 


X 


10 


-4 


7.26 x 10+ 8 


26.4 


0.891 


1.02 


1 


.42 



Note. — The optical depth at the turnover, r m , is calculated for a homogeneous sphere. The factor by which the extrapolated 
optically thin flux at the turnover, S m , over-predicts the observed peak flux, So, is also calculated for the geometry of a 
homogeneous sphere. 
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